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For deuteron electromagnetic form factors, perturbative QCD(pQCD) pre- 



(-H , diets that Gq becomes the dominate heficity amplitude and that G^ and 

G+_ are suppressed by factors Aqcd/Q and Aq CD /Q; 2 at large Q 2 , respectively. 
We try to discuss the higher order corrections beyond the pQCD asymptotic 
predictions by interpolating an analytical form to the intermediate energy 
region. From fitting the data, our results show that the helicity-zero to zero 
matrix element Gq dominates the gross structure function A(Q 2 ) in both 
of the large and intermediate energy regions; it is a good approximation for 
to ignore the higher order contributions and the higher order correc- 
tions to G^q should be taken into account due to sizeable contributions in the 
intermediate energy region. 
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1 Introduction 



It was found a long time ago that the traditional meson-nucleon picture can 
not explain the form factors of the deuteron as the momentum transfer Q 2 > 1 
GeV 2 W. It means that the fundamental degrees of freedom of QCD, the quark 
and gluon degrees of freedom, must be taken into account. However, A pure per- 
turbative QCD (PQCD) calculation' 2 ^ shows that the theoretical prediction is much 
smaller than the data at currently accessible energies, although it may be correct 
in very large Q 2 . To explain the deuteron form factors in the intermediate energy 
region, we have suggested a QCD-inspired model for the helicity-zero to zero matrix 
elements Gq in the light-cone framed, which should be the dominant amplitude 
from PQCD predictions' 4 '. This model can explain the data of the deuteron electro- 
magnetic structure function A(Q 2 ) and shows that Gq is already dominant at Q 2 
of 1 GeV 2 . Furthermore, it was found that G| can not be neglected in the form 
factor G M (Q 2 )^'^ and G\_ plays an important role in Gq(Q 2 )^. Neglecting them 
will result in contradictions with both the data and the conventional meson-nucleon 
picture in the low energy region. 

Perturbative QCD (PQCD) predicts that Gq becomes the dominant helicity 
amplitude at large Q 2 and that G+ and G+_ are suppressed by factors Aq C d/Q 
and Aq CD /Q 2 , respectively, neglecting G\ and G\_ contributions at large Q 2 , we 
have the relation approximately, 

G c : G M : G Q = (1 - Hij) : 2 : -1, (1) 

where rj = Q 2 /AM 2 and M is the mass of the deuteron. However, the helicity- 
flip amplitude G+ and G+_ contribute to Gm and Gq because of the kinematic 
enhancement in the intermediate energy region. In order to explore the role of G\ 
and G+o, we have tried to expand them to the second order in Aqcd/A^, according 
to the QCD predictions at large Q 2 . This expansion can connects smoothly PQCD 
predictions in the high energy region with traditional nuclear physics predictions in 
the low energy region. It was shown that the second order contribution strongly 
affects the behavior of Gq in the intermediate energy region. At large Q 2 , the ratio 
of form factors (1) is slightly modified. 

Following this approach, G+ and will be expanded to higher orders (beyond 
the second order) in Aqqb/M according to the PQCD prediction at large Q 2 . In 
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order to explore the role of higher order contributions we discuss the possibility to 
interpolate an expression for G+ and G|_ to the intermediate energy region in this 
paper. It is worthwhile to unify the predictions for the deuteron form factors from 
the low energy to large energy region. 

A general consideration based on the PQCD predictions for the deuteron from 
factors is given in Sec. 2. As an example, a phenomeno logical analytic form including 
the higher order corrections in Aq C d/^ is suggested in Sec. 3. The numerical results 
and summary are presented in Sec. 4 and Sec. 5, respectively. 



2 A General Consideration Based on the PQCD 
Predictions 

For the deuteron case, the matrix element of the electromagnetic current is 
defined as 

G% x = (P'X' | r | PA) , (2) 

where Q 2 = —{P' — P) 2 and \PX) is an eigenstate of the deuteron with momentum 
P and helicity A. In the standard light-cone frame (LCF), defined by Ref.[7] q + = 
0, q y = 0, and q x = Q, the charge, magnetic, and quadrapole form factors can be 
obtained from the plus component of three helicity matrix elements: 
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In terms of Gq, Gm, and Gq, the Rosenbluth cross section and the tensor polariza- 
tion T 20 for elastic ed scattering can be expressed as 
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where A(Q 2 ) and B(Q 2 ) are given by 



and 



A(Q 2 ) = Gl + \r)G 2 M + l V 2 Gl (6) 



Bm = ^v(l + V)G 2 M . (7) 



It was shown M that, in LCF, the helicity-zero to zero matrix element Gq would 
be the dominant helicity amplitude at large Q 2 for elastic ed scattering from the 
PQCD predictions. It means the Gq dominance in the structure function A(Q 2 ). 
It was also argued' 8 ! that the dominance of Gq Q begins at Q 2 ^> 2MAqcd ~ 0.8 
GeV 2 but not rj » 1. Thus 2MAqcd is a scale of validity of PQCD predictions 
and the quark and gluon degrees of freedom in the deuteron should be taken into 
account to solve the problem that the experimental results of A(Q 2 ) are in sharp 
disagreement with the meson exchange calculations for Q 2 > 0.8 GeV 2 ' 1 !. To make 
detailed prediction for Gq , we have suggested a QCD-inspired model' 3 ! in the region 
of Q 2 > 1 GeV 2 based on the reduced form factor method' 9 ! , which fit the data well. 

PQCD predicts that G+ and G\_ are suppressed by factors Aq C d/Q and Aq CD /<5 2 , 
respectively. However, in the intermediate energy region, Gq dominates the charge 
form factor Gc, but not Gm and Gq. As shown in Eq.(3), while rj < |, the G\ 
contribution to Gm and Gq are enhanced by a factor and G+_ contribution to 
Gq is enhanced by a factor Although G+ and G|_ are suppressed for dynamic 
reason, they contribute significantly to Gm and Gq because of the kinematic en- 
hancement. Without these contributions, the predicted form factors, except for Gc, 
are in sharp disagreement with the data. In order to explore the role of G+ and 
G\_, we interpolate a general expression based on perturbative QCD predictions, 

Gto = ^9 + o(v) G+ (8a) 



G%- = ^j9+-(v) Gto , (8b) 



where g+o(v) and g + -(rj) are any functions of r] with r] = -^p- Obviously G+ and 
are suppressed by factors Aq C d/Q and Aq CD /Q 2 as long as g+o(r)) and g + _(rj) 
satisfy the following condition, 
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9+o(v), 9+-(v)^0(l) as r]^oo . (9) 
Substituting Eqs.(8) into Eqs.(3), one can get 



°< = 2p+ (2Ui) l(1 - !"» + 1 *»<"> + 5< 2 " - l) T n 9 - ( " )|G » • (10o) 



1 2n — 1 1 

G " = ^TT)I 2 + V 9+o( " ) " 5? J+ - ( " )|G » 

and 

1 r 1 1 / n V+l 



G c 



2^TT) I" 1 + ~ v ^ ~ V 5+ - (?7)]Go+0 • (1 ° c) 



Now we discuss the constriants on g+o(v) an d g+- (v) from the experimental data. 
As we know, G M (Q 2 ) changes sign at Q 2 = Ql ~ 2GeV 2 ^ (or rj = r) = ^ ^ 
0.13). The dominance of Gq can not explain this point (see Eq.(3b)) and at least 
the second term in Eq.(3b) should be in the same order of the first term to cancel it 
in order to fit data of Gm- That means g+oiv) is nonzero. If we first keep g+-(rj) = 
in Eqs.(10), then g +0 (r] ) can be determined by the zero at Gm(t]o)> which turns out 
to be g+o(r]o) = jz^- m this case, Gq is negative at rj = rj where PQCD begins 
to be valid. Thus there must be a node in the region Q 2 < lGeV 2 since Gq is 
positive at the origin experimentally. The theoretical prediction is contrary to the 
experimental data without G+_ contribution. Therefore the predicted form factors 
are in sharp disagreement with the data without G\ Q and G\_ contributions and the 
existence of non-zero g+o(i]) and g+-{rj) is necessary. In addition to the constrant 
(9), g+o(i]) and g+-{rj) should astisfy 



/ x 2?7o — ^+-(770) 
9+o(Vo) = ; 7T I 11 ) 

1 — Z.TjQ 



to ensure Gm(vo) — and 



2 V g +0 (r ] )-2r ] 2 <( V + 1) g + _(rj) (12) 

to keep Gq positive at any momentum transfer. In particular, as rj = rj Eqs.(ll) 
and (12) give the constraints on g+o(vo) an d 9+-(.Vo), 

9 + -(Vo) > (13) 
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and 

ff+ofa)) < T% (14) 
Eqs.(9), (11) and (12) are three constraints on the functions g+o(v) an d 9+-(v)- 

3 A Phenomenological Example 

As mentioned in the section 2, g + -{rj) ^ is important to keep Gq > although 

A 2 

G +- = ^9+-(v)Go is suppressed by the higher order factor A 2 QCD /Q 2 (= -f^- ■ i). 
However the G+ is the first order correction which makes the zero in Gm at ~ 
1.85GeV 2 . We have expanded G\ Q and G\_ to the second order in Aqqd/M in 
Ref. [6] and numerical results show that the second order contribution to G+ plays 
an important role in the intermediate energy region. In this paper we introduce an 
exponential form phenomenologically as an example, 

9+o{v) = f exp( — ) (15a) 

V 2 ^ 



and 

9+ -(v) = fexp(-^=) (156) 



to interpolate the higher order corrections. Obviously, the exponential form (15) 
satisfies Eq.(9) and is consistent with the perturbative QCD prediction at the large 
transfer momentum region. Thus Eq.(15) is enable us to make an analytical evalu- 
ation to the deuteron form factors. 



4 Numerical Results 

We input the parameters b and c, and determine / by the zero in Gm- For a 
certain c, the obtained Gq increases proportionally with 6. We can fix b to connect 
our predictions with the data smoothly. To retain good convergence, we constrain 
o^7= and being smaller than unity. On the other hand, we restrict b and c 
to be positive to keep the exponential damping as Q 2 goes to infinity and it is a 
resonable assumption, after taking into account the higher order corrections. For 
c = 0.0, 0.5 and -0.5, the predicted B{Q 2 ), G Q (Q 2 ) and T 20 {Q 2 ) are shown in figs. (1- 
3). Experimental data are taken from Refs. [1,10-12]. To smoothly connect with the 
data of Gq, the parameter b should be 1.1, 1.3, and 0.8, / = 0.37,0.51, and 0.30, 
respectively. As argued above, c = —0.5 should be abandaned. While c = 0.5, 
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the corresponding / is 0.51, which is too large to keep smaller than unity as 
Q 2 > 1. The parameter c = 0.0, b = 1.1, / = 0.37 is an appropriate choice. 

5 Summary 

Based on perturbative QCD predictions at large momentum transfers we have 
tried to discuss the corrections to deuteron form factors beyond the second order 
in Aqcd/M.A general consideration is given by introducing functions g+o(v) an d 
g+-(r)) and the data at the present energy region put three constrants on the func- 
tions g+o{v) an d 9+-(v)- I 11 order to explore the role of higher order contributions we 
suggest an exponential form for g+o{v) an d g+-(v) as an example. We conclude that 
(1) the helicity-zero to zero matrix element Gq dominates the gross structure func- 
tion A(Q 2 ) in both of the large and intermediate energy regions. A QCD-inspired 
model can describe this matrix element well®. (2) and G+_ contributions are 
important in determining form factors Gm and Gq; (3) By fitting the data, we get a 
set of parameters in g+o(v) an d 9+-(v), c = 0-0? b = 1.1, / = 0.37, which can describe 
Gm, Gq and T20 appropriately, the parameters c = 0.0 indicates that it is good for 
G\_ in the intermediate energy region to take the asymptotic behavior which were 
predicted by pQCD, and the higher order (beyond the second order) contributions 
are negligible; (4) The higher order corrections to G+ should be taken into account 
and they make sizeable contributions in the intermediate energy region. 
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Figure Captions 



Fig.l. Structure function B(Q 2 ). The dashed dotted line corresponds to the Paris 
potential calculation. Experimental data are taken from Ref. [10]. 

Fig. 2. The form factor Gq. Experimental data are taken from Ref. [11]. 

Fig. 3 The tensor polarization T 20 with scattering angle 9 = 70°. The dashed dotted 
line corresponds to the calculation with Paris potential. Experimental data 
are taken from Ref. [11] 
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